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VECTOR ANALOGUES OF THE MAGGI-RUBINOWICZ 
THEORY OF EDGE DIFFRACTION 

Robert Meneghini, Peter Shu and John Bay 
ABSTRACT 

The Maggi-Rubinowicz technique for scalar and electromagnetic fields can be interpreted as a 
transformation of an integral over an open surface to a line integral around its rim. Using this trans- 
formation, Maggi-Rubinowicz analogues are found for several vector physical optics representations. 
For diffraction from a circular aperture, a numerical comparison between these formulations shows 
the two methods are in good agreement. To circumvent certain convergence difficulties in the Maggi- 
Rubinowicz integrals that occur as the observer approaches the shadow boundary, a variable mesh 
integration is used. For the examples considered, where the ratio of the aperture diameter to wave- 
length is about ten, the Maggi-Rubinowicz formulation yields an 8 to 10 fold decrease in computa- 
tion time relative to the physical optics formulation. 


VECTOR ANALOGUES OF THE MAGGI-RUBINOWICZ 
THEORY OF EDGE DIFFRACTION 

Introduction 

The theory of the edge diffraction wave originated with Thoma: Young, who observed that 
light diffracted through an aperture could be interpreted as an unperturbed incident wave com* 
bined with a wave disturbance arising at the aperture rim [1, 2] . Starting with the Helmholtz 
representation, the idea was formulated in a mathematically rigorous way by Maggi and later by 
Rubinwicz [3] . A comprehensive review article by Rubinowicz ( 1 ] traces the development of the 
theory up throu^ the work of Miyamoto and Wolf. 

Although the approaches of Maggi and Rubinowicz are different [1,2,4], both lead to the 
I reduction of a surface integral to a line integral (usually corresponding to the physical edge of the 

scatterer) plus geometrical optics terms. Miyamoto and Wolf [5] were the first to demonstrate 
that such a transformation could be carried out exactly for an arbitrary incident field. Although 
their work centered on the transformation of the Kirclmoff integral, it was later shown that a sim* 
ilar procedure could be carried out for the Rayleigh and Sommerfeld integral representations [6] . 
In a more recent paper, expressions have been derived for the field behavior near the geometrical 
optics shadow boundary [7]. Apart from the importance of the results themselves, they provide 
information on the relationship of the Maggi-Rubinowicz technique to the geometrical theory of 
diffraction (GTD) and to the Braunbek approximation [8, 9]. 

M(»t of the papers cited above deal with the scalar formulation. For electromagnetic diffrac- 
tion problems, vector analogues of the Maggi-Rubinowicz technique have been investigated by 
several authors [1,10, 11]. The application of these ideas to antenna problems, however, does 
not seem to have attracted much attention and most of the practical results [12-15] have been 
i presented for typical applications in optics. An exception is a paper by Gordon [ 16] , who ob- 

I tained scattered fields from a planar reflector. In this work, however, only scalar fields were 

1 considered. 
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In this paper we obtain the Maggi-Rubinowicz (M-R) analogues to three commonly used 
physical optics (P. O.) formulations. Two of these have been derived previously by somewhat 
different methods than the one given here [1,11]. We next present some numerical examples 
comparing one of the M-R formulations with the corresponding P. O. representation. For the 
simple examples considered, we show that a variable mesh integration can be used to circumvent 
errors that normally occur when the observer approaches the geometrical optics (g. o.) shadow 
boundary. 


Some Maggi-Rubinowicz Representations 

Using an time convention, we list below three P. O. representations for the scattered 

field. 


Ep = \7x f (n'xEjGdS'- — (77 - + k2) | (n'xHJGdS' 

I 


E, = 2 7x 
^1 


.) G dS' 


I 2 »- -r— (77 - + k2) 


Koe 




)GdS' 


( 1 ) 


( 2 ) 


(3) 


where G is the free space Green’s function 

-ikix-x'i 

G= ^ ~ ~ 

4ir|x - x'l 

and where x, x' denote the coordinates of the observation point and source point respectively. For 
all three formulations, the surface of integration, S, is open. Denoting a closed surface as the 
union of the surfaces S and S, where the respective integrands are assumed to be zero on S, we 
choose n' to be the inward normal to S^. For the fields appearing in the integrands we have used a 
subscript ‘a’ to denote the approximate nature of the fields. To simplify the notation we will omit 
the subsaipts in the subsequent equations. 

Equations (2) and (3) apply to the problem of diffraction throu^ an aperture in a perfectly 
conducting plane screen [ 17] . The observer is limited to the shadow half space, i.e., the half space 
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excluding the sources of the incident field. The relevant approximations are E, ■ H, • where 
Eq, are the incident fields with the surface of integration being taken over the aperture. The com- 
plimentary problem of the screen can be found directly from E| , JI 2 Babinet’s principle [17] . 

Equation (1) is often referred to as the Franz formulation [ 18, 19]. It is identical to the ex- 
pressions given by Stratton [20] and Kottler [21, 22] despite differences in appearance from these 
last two. For scattering from a perfectly conducting (and in general non-planar) object, n' x E 0 
and n’ X 2(n' x H^) where the open surface S corresponds to the illuminated portion of die 
scatterer. 

To find the M-R analog to Ep we first relate Ep to the vector Kirchhoff representation, Ej^, 
where [23, 24] 


Sk*J (E(a'*7'G)-G(8'-7')E|dS' 


(4) 


the integrand can be rewritten as 


2E(8'-V'G)-(n'*V')GE 
Using v' • E s 0 and the vector identity 

(b • 7)a “ 7(a • b)-(a ‘Vjb - a x(Vxb)-bx(7xa) 


then 


(n' • V')(GE)«-G(a'x(7'xE))-ft'x(7'GxE) + ft'(E • 7'G)»X *8' 
The quantity ^ is a dyad, the ji component of which is 

where 5jj is the Kronecker delta function. 

Using the relationships 


and the identity [ 1 ] 


V'xE »ici>MH 

n'x(7'GxE)»(8' E) 7'G-(8' *7'G)E 


I I - 8'dS'» |(GEx«')d£' 


QP roOR 
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where fi' is a unit vector along the closed contour C directed so that fi' and n' are related via the right 
hand rule, then (4) becomes 

EK = f ((n'xE)xV'G + iwMG(n'xH) + (n'-E)7'Gl dS' (5) 

Js 

+ ^ G(ExK')dfi' 

The surface integral in this equation is often referred to as the Stratton*Chu integral, which we de- 
note by Eg_c- It should be noted that the derivation of Ej_^. from Ej^ is sometimes done under the 
assumption of a closed surface, e.g. [ 24 J , so that E * Ej^. The rim integration in (5) is therefore 
a consequence of assuming the surface to be open. 

The relationship between Eg_j. and Ep is given by [ 17] 

= V(tGH-8'd«' (6) 

This equation has been derived by Stratton [20] and Kottler [21, 22] and is obviously equal to the 
Franz representation, (1 ). 

In order to express Ep solely in terms of line integrals we first write equation (4) in the form 

where the implicit sum on i runs from 1 to 3 and where 

Vi(x') = Ei(x') 7'G-G7'E:(x') 

Following Miyamoto and Wolf [5] we notice that since 7' • Vj = 0, a vector potential Wj can 
be defined such that V. = 7'xW., Assuming that W. has only a finite number of singularities that are 
excluded from the region of integration, then Stokes theorem yields 

" *1 {£ Si • '■ ‘“’j *''* 

where the notation Ej^ is used to distinguish the P. O. formulation of (4) from the corresponding 
M-R formulation. 
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The second term of (7) represents integrations around the singularities of W|. The unit vector 
fi* around the contour Cj is directed in the opposite sense of around the rim c. The vector poten- 
tial W| is given by [S] 

W,(x,5') = -G r dM + W^ 


where, as before, x,x' are the vectors from die origin to the observer and source point respectively 
and r - (x - x')/| X - X Miyamoto and Wolf have shown that is zero if the field satisfies the 
Sommerfeld radiation condition. These authors have shown, moreover, that if each component of 
the electric field can be written in the form Ej « A^e^, dien W| can be represented as an asymptotic 
power series in the wavenumber k, 

A A 

Wj = -GEi +0(k-») 

» 1 - r* p 

with 


p= 


For plane or spherical incidence waves all but the first term of the expansion for W| vanish. 
Writing (7) in dyadic notation, with W = Xj Wj, 

Ik’ ^ W *fi'd£' + ^^ W •K'd£' 

then from (S) and (6), the M-R analogue of the Franz representation, Ep, becomes, 

Ep= W £'-GEx 2'+— VG(H -2') d£'+y i W-£'dk' 

^ i«/i ~ j 


( 8 ) 


A similar formula has been derived by Rubinowicz [ 1 ] . 

The M-R analogues of (2) and (3) can be obtained in a straightforward manner by using a method 
d^cribed by Jones [17]. For an aperture in a perfectly conducting plane screen in the z«0 plane, we 
can write the scattered field in the diffraction half space by means of ( 1), Le., Ep(x), with x = (x,y,z). 
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For an observer at the image point, we employ the same formula obtaining Ep(x f). x | (x,y,-z) 
where the Greens function in the integrand now becomes 


gikf, 


G(x,.x')»- 
"• 4»rr| 

where rj « |x j - x' | is the distance from the source to the image point. Since the observer at the 
image point is outside the original volume of integration then Ep(x |) » 0. Since the source points 
are located in the z = 0 plane, then G(x, x')* G(X|, x), x • ^G(x, x') = x • VG(Xj, x'), y • VG(x, x')= 
y '^GfXj.x'l.z • ^7G(x, x')“-z • VG(X|,x'). Using these facts it is not difficult to show ( 171 that 
by adding the z components of Ep(x) and ]|p(X|) and subtracting the x, y components of Ep(X|) 
from Ep(x) we obtain equation (2). Equation (3) is obtained by adding the x, y components and 
subtracting the z components. We can express these operations throu^ the combinations 

I -EpCxl-Q - 2zz) *Ep(Xj) (9) 

i •Ep(x)+(|:-2zz) -Ep(x,) (10) 

A A . A /X A A 

i * XX + yy + zz 

where (9) yields (2) and (10) yields (3). Now by replacing Ep by Ep in the above equations (where 
Ep is the M'R analogue of Ep) we directly obtain the M-R analogues of Ej , E 2 that we denote by 
E' E' • 

£p£2- 

E'|- ji |d, •l'-2zG(Ext')-z|(ir*^^ D| 8'df (ll» 

Ej- £|Dj.8’ + 2G(I-22).(Ex«'),J-70(H.8)|,^|^ Dj -S' dt' 


where 

D, =I •W-(I-2zz) -W' 

s 1 s s s s 

D 2 = I *W + (1-2zz)W' (14) 
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with 

W-Jf(X,X')-S,W,(X.K') 

W-.W(x„x')-8,W,(x„xO 

Equation (I I) has been obtained by Kaiczewski ( 1 1 ) by generalizing the work by Mardiand and 
Wolf [6] to the vector case. 


Results 

To compare the M*R and die P. 0. formulations we numerically compute the field diffracted 
from a circular aperture. To do this we use (2) for the P. O. and (1 1 ) for the correqxmding M*R 
formulation. For plane and spherical wave incidence (the two cases that will be considered) the 
vector potentials simplify to 

Assuming the perfectly conducting plane screen to lie in the z * 0 plane with the sources of the 
incident field located in the z < 0 half space, then for observation points z > 0, (1 1 ) can be written 

E', - f g|e,(A, •«') + zE,(A2-h-2i(Exfi').z|d£' + g.o. (15) 


where 






2 


A A 


rjxp 


1-^ 


rxp 

; — X — jr 

1 -r *p 


(16) 


The term labeled g. o., which equals the second term of ( 1 1 ). is simply the vector geometrical optics 
field. This is determined by tracing the incident rays through the aperture, keeping track of the phase 
and amplitude, and preserving the vector nature of the inddent field. That this term is a g. o. field 
follows from the work of Miyamoto and Wolf (51 and the fact that W' has no singularities over the 
aperture since 1 - T| * p is never zero there (6) . 
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As the observation point approaches the g. o. shadow boundary the singularity of W approM;hes 
die aperture tim. When the observer is direcdy on the shadow boundary a separation between the g.o. 
term and the rim integration is no longer possible. However, for the numerical computation we can 
use the fact that even with a high sampling rate in the observation space, it is highly improbable that 
the observer will lie so close to the shadow boundary that the value of the integrand will result in 
computer overflow. To properly account for the contribution to the inte^l near this singular point, 
the sampling rate along the rim must be increased. To accomplish this a variable sampling grid is used 
for the integration along the aperture rim. The sampling rate is determined so that the integrand 
clianges by some smaU fraction of its value in moving from one point to the next. For the observer 
near the shadow boundary this procedure yields a fine sampling grid for that section of the aperture 
edge near the singularity and a more coarsely sampled grid along the remaining portion ot the edge. 

It is evident that this procedure circumvents rather than solves the underlying analytic problem. 
To obtain uniformly valid expressions as the observer approaches the shadow boundary it would be 
necessary to generalize the work of Otis, et al (7) to the vector case and to arbitrary aperture shape*. 

In figures I -5 are shown comparisons between the M-R (solid line) and the P. O. formulation 
(X). We have chosen throughout a frequency of 4 GHz and an aperture radius of 0.4 m. The screen 
and aperture are located in the z ■ 0 plane. The diffracted field has been computed for points in the 
X > 0, z > 0 quadrant of the x-z plane. The absissa ot the graphs corresponds to the angle, in radians, 
measured from the z axis. The ordinate corresponds to the magnitude of a particular component of 
the field which is given in units of volts/meter or millivolts/metcr. For the M-R results, a variable 
mesh integration around the aperture rim was used to generate all but the final figure. 

For an incident monochromatic (dane wave of the form 

Ej-ec**i 

whew 

e ■ ((k • z) X -(k • x) z) 

A ^ 

k * z ■ cos(ir/8) 

A ^ 

k • x«sin(ir/8) 
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we plot in flgures I and 2 the magnitude of die E, components of the Held. The distance, r^, 
from the center of the aperture to the observer is S m. It can be seen that the results from the P. O. 
and M*R representations are in good agreement. 

In figures 3 through S we assume an incident plane wave of the form 

E‘-0'f(d')e“‘Vr 


where r is the distance measured frrnn the source poii.t. Hie distance from the source to the center 
of the aperture is taken to be 3 m. The polar angle 6 ' is measured with respect to the z' axis where 
X • z’ * cos(ir/8), y z' ■ 0 z • 2' ■ -sin(ir/8) and where the unit vectors x, y correspond to the x, y 
axes located in the plane of the aperture. Hie unit vector associated with 6 * is denoted by 6 '. For 
the numerical computations we have assumed an isotropic source so that f(9') ■ 1 . 

In figure 3, the magnitude of E, is plotted for a far field observer, r^ > 5 m. A fixed sampling 
grid was used to generate figure S. Hie two discontinuities that occur near die opposite points of 
the shadow boundary in this figure are eliminated by using the variable mesh integradon (figure 4). 

The m^or advantage of the M-R relative to the P. O. representation is the smaller amount of 
computer time needed to calculate the diffracted field. For the parameters chosen here, with a di- 
ameter to wavelength ratio, D/X, of about 10, the time needed for the M-R calculation is about a 
factor of 8 to 1 0 less than that required for the P. O. calculation. This factor in savings would in- 
crease for larger DA. Foi more complicated geometries, e.g., non-planar reflectors, the computa- 
tional savings is not expected to be as great since the calculation of the g. o. terms would be more 
time consuming than in the simple aperture problems considered here. 

It should be mentioned that for simple aperture shapes such as circular and rectangular, the 
Frauenhofer approximation can be used to obtain closed form expression for the M-R and P. O. 
representations ( 1 5) and (2). In particular, the results for a plane wave incident upon a circular aper- 
ture have been shown to be in good agreement for i>hsi. rvation points outside the geometrically illum- 
inated region (Appendix A). The M-R solution, however, is discontinuous near the shadow boundary. 
For this case the variable mesh integration is not a^ plicabie since the integration is performed prior 
to any numerical computation. 
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Applying the techniques used by Keller et al. ( 9] and Miyamoto and Wolf [ S 1 to equations ( 
and (I S), it can be shown that the flrst terms from the asymptotic expansions in k arc in agrecmcut. 
Hiis result holds for arbitrary aperture shapes and arbitrary Helds of incidence. The expansions that 
were used, however, are not valid near the shadow boundary or near a caust.. 
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Figure 1 . M-R (solid line) and P. O. (x) solutions for |E^| , r^ » Sm; 
plane wave incidence. 
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Figure 2. M-R (solid line) and P. O. (x) solutions for lE^I , Fq = 5m; 
plane wave incidence. 
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Figure 4. M-R (solid line) and P. O. (x) solutions for lE^I , rj, = 5m; 
spherical wave incidence. 
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spherical wave incidence. Uniform sampling mesh along the aperture ^ 

rim chosen for M-R solution. ; 
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APPENDIX A 


Compuison of P. O. and M*R in the Frauenhofer Region 

For simple aperture shapes, closed form expressions can be obtained for the P. O. and M-R for- 
mulations in the Frauenhofer zone. In particular, we compute E ] > Ej and E|^ for the case of a plane 
wave inddent upon a drcular aperture. We place the screen in the z « 0 plane with the sources of 
the incident field in the z < 0 and the observer in the z > 0 half space. The inddent electric field is 
assumed to be. 


Ej»Se‘k-i' 


where 


To compute E| from the formula 


e = [(k • z)x-(k • x)zl 
k • z = cos(ir/8) 
k • x = sin(jr/8) 
k = kic 

Jl, =2 Vx|^(7xE)GdS' 


(Al) 


we use the vector relations 

VxI(zxE)G] = ^7Gx(zxE) + G7x(zxE) 

VGx(zxE) = (7G • E)z-(VG • z)E 

Noticing that 7x(zxE) ° 0 since E is a function only of the primed coordinates then 

E, »2| |z(7G • E)-(VG • z)e| dS' 

Furthermore, 

VG = f(ik-I/r)G~ikfoG 

e^'^o .. , 

G ~ I 
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Substituting this result into (A1 ) and (A2), then 


ilc r 

iae*^*^o A 


Ej (k • z)J,(kaf)(roXy) 

$ ’’o 

£ 2 ° (k-^ro) • zJ,(kar)e 

where a is the radius of the aperture, and J| is the Bessel function of order 1 . 

To find the M-R analogue of E | we begin with equation (1 5) of the text, 

E'j » ^ g|Ej(A, • fi') + zEj(A 2 • 8')-2z(Ex«') • z|d«' + g.o. 


where 


E,»iE,+9Ey 

A A 

r,xp 


A A 

A A 

rxp 

A s 4 ■ ■ _ I... ■ ■ 

1*«A A *AA 

1-r, -p 1-r p 


We again make use of the Kirchhoff assumption 

Et(x') = x(k • z)e‘!i'i’ 


E_(x')*-(k • x) e‘l5 'i 

^ C *s» 


uiu the relations 


(A3) 


(A4) 


^ t 

p«k 

S' = (-sin 0', cos 0 ', 0) 

(fx p) • 8' ~ -sin 0' l(k • z) sin 0 cos 0 - (k • y) cos 9) 

+ cos 0'l(k • x) cos 0 - (k • z) sin 0 cos 0) 
r • p (k • x) sin 0 cos 0 + (k • z) cos 6 

where the observation point (x, y, z) has been expressed in terms of a spherical coordinate system 
so that xscos 0 sin 0, y ® sin 0 sin d, z = cos 0. The expressions for (r|Xp) • I f| • p follow from 
the above formulas by replacing cos 0 with -cos 9. 

After the resulting integrand is simplified we encounter the integral, 

A ^ 

gikal(k.x)co8«. -cos(<.-o')8ine, ^ cos(0 - 0') + a . cos 0') d0' 

I • ^ 
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Using the transformation -f u, described above, the integral can be shown to reduce to 

( 0 | cos (0 • u) «2 cos u) 2iri J| (ka{) 
where Oj , Oj are constants and di u are defined above. 

After rearranging terms, E | can be expressed by the complicated formula. 


I ra A 

ll ■ —7- J|0cef) (-X d, ♦ i(fi2 - (k • i) a)l ♦ g. o. 

»o » 


(AS) 


where 


(fc-fo) * 


1 


dl ■ “ (k • z) (fo • z) j(fo • z)2 Ur, • y)* - o(fo • x)l (k • x) a a} 


^2 “ ” * *)* (1^ * *)(^ • x)o-o(k • z) (o(r, • x)-(r, • y)*l| 

o-l-(k-xHfo-x) 

7«(1 -(k • x)(r, • x))2-(f, • z)2 (k • z)2 

In figures Al, A2 we have plotted the magnitude of Uie E, components of the field for both 
the M*R (solid line) and the P. O. (x) formulations. The observer is taken to be at a distance of 10 km 
from the center of die aperture. We assume, in addition, an incident plane wave identical to that con* 
sidered in the text. All other parameters are the same as those ^en in the text: f * 4 GHz and an 
aperture radius equal to 0.4 m. The geometrical optics term of (AS) lias not been included in the 
numerical computation. In consequence, the M*R formulation, as plotted, is incorrect in the small 
angular region at the cent^ of the main lobe. 

However, even if the g. o. term were added to the solution, the discontinuity at the shadow 
boundary would persist. It therefore appears necessary to use a numerical integration with the vui* 
able grid method for those points within the g. o. illuminated region. For obsnvation points outside 
this region it can be seen from the figures that the two solutions yield nearly the same results. 
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Figure Al. M>R (solid Une)and P. O. (x)solutioiis for lE^I , fg • 10 km: 
plane wave incidence; Frauenhofer solution. 
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Figiore A2. M*R (solid line) and P. O. (x) lolutioiis for |E,| . ■ 10 km; 

plane wave incidoioe; Frauenhofer solution. 
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